Kubo formula for Floquet states and photoconductivity oscillations in a 2D electron 

gas. 



Manuel Torres 1 ^ and Alejandro Kunold 2 -0 

1 Institute/ de Fisica, Universidad National Autonoma de Mexico, 
ly-j | Apartado Postal 20-364, Mexico Distrito Federal 01000, Mexico 

. 2 Departamento de Ciencias Bdsicas, Universidad Autonoma Metropolitana-Azcapotzalco, 

Q ■ Av. San Pablo 180, Mexico D. F. 02200, Mexico 

(SJ ; (Dated: February 2, 2008) 

The recent discovery of the microwave induced vanishing resistance states in a two dimensional 
electron system (2DES) is an unexpected and surprising phenomena. In these experiments the 
magnetoresistance of a high mobility 2DES under the influence of microwave radiation of frequency 
C^") ' u-> at moderate values of the magnetic field, exhibits strong oscillations with zero-resistance states 

(ZRS) governed by the ratio lu/uj c , where co c is the cyclotron frequency. In this work we present 
a model for the photoconductivity of a two dimensional electron system (2DES) subjected to a 
magnetic field. The model includes the microwave and Landau contributions in a non-perturbative 
exact way, impurity scattering effects are treated perturbatively. In our model, the Landau-Floquet 
states act coherently with respect to the oscillating field of the impurities, that in turn induces 
transitions between these levels. Based on this formalism, we provide a Kubo-like formula that 
takes into account the oscillatory Floquet structure of the problem. We study the effects of both 
short-range and long-range disorder on the photoconductivity. Our calculation yields a magnetore- 
sistance oscillatory behavior with the correct period and phase. It is found that, in agreement with 
experiment, negative dissipation can only be induced in very high mobility samples. We analyze the 
dependence of the results on the microwave power and polarization. For high-intensity radiation 
multi- photon processes take place predicting new negative- resistance states centered at u>/u) c = 1/2, 
and lu/lo c — 3/2. 

O ' I. INTRODUCTION. 

l cj; 

Two-dimensional electron systems (2DES) in a perpendicular strong magnetic field have been extensively studied 
in relation with the quantum Hall effect. Recently, two experimental groups^^*^, reported the observation of a novel 
phenomenon: the existence of zero-resistance states (ZRS) in an ultraclean GaAs / Al x Ga\- x As sample subjected 
to microwave radiation and moderate magnetic fields. The magnetoresistance exhibits strong oscillations with ZRS 
\ governed by the ratio e = lo/uJc, where lo c is the cyclotron frequency. According to Zudov etal., the oscillation 
f**- ■ amplitudes reach maxima at e — lo/lj c = j and minima at e — j + 1/2, for j an integer. On the other hand Mani 
etal. reported also a periodic oscillatory behavior, but with maxima at e = j — 1/4 and minima at u/u) c = j + 1/4. 
"^t" , Additionally experimental work appeared recently^SiLSiS. 

In spite of a large number of theoretical WO rksi2iiA4i2ii2ii2ii£ii£iiLi2ii2iSi ; a complete understanding of the new effects 
in 2DES induced by microwaves has not yet been achieved. A pioneering work put forward by Ryzhii^iii predicted 
the existence of negative- resistance states (NRS). Durst and collaborators^ found also NRS in a a diagrammatic 
calculation of the photoexcited electron scattered by a disorder potential. A possible connection between the calculated 
*^ • NRS and the observed vanishing resistance was put forward in reference^, noting that a general analysis of Maxwell 
equations shows that negative resistance induces an instability that drives the system into a ZRS. Whereas some of 
the modelsi2*ii*i2iil are based on an impurity assisted mechanisms, there are alternative explanations£*i2i2S in which 
the leading contribution arises from the modifications of the electron distribution function induced by the microwave 
radiation. These models, as well as some of their predictions remain to be tested experimentally. 

In this work we present a model which includes the Landau and radiation contribution (in the long-wavelength 
limit) in a non-perturbative exact way. Impurity scattering effects are treated perturbatively. With respect to the 
Landau-Floquet states, the impurities act as a coherent oscillating field which induces the transitions that proved to 
be essential in order to reproduce the observed oscillatory behavior of the magnetoresistance. Based on this formalism 
a Kubo-like expression for the conductance is provided. Our results display a strong oscillatory behavior for p xx with 
NRS. It is found that p xx vanishes at e = u>/u> c = j for j integer. The oscillations follow a pattern with minima at 
e = j + 5, and maxima at e = j — 5, adjusted with S » 1/5. The model is used to test chirality effects induced by the 
magnetic field, calculations are carried out for various .E-field polarization's. Finally, we explore the nonlinear regime 
in which multi-photon processes play an essential role. 

The paper is organized as follows. In the next section we present the model and the method that allow us to obtain 
the exact solution of the Landau-microwave system, as well as the perturbative corrections induced by the impurity 
potential. In section UTTl we develop the formulation of dc electrical linear response theory valid in arbitrary magnetic 
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field and microwave radiation. A discussion of relevant numerical calculations is presented in section (|IV|I . The last 
section contains a summary of our main results. Details of the calculations are summarized in the appendices. 

II. THE MODEL. 

We consider the motion of an electron in two dimensions subject to a uniform magnetic field B perpendicular to 
the plane and driven by microwave radiation. In the long-wave limit the dynamics is governed by the Schrodinger 
equation 



<9<f 

ih— = iJ* = 
dt 



H {BiU} + V(r) *, (1) 



here H^b.u} is the Landau hamiltonian coupled to the radiation field (with A — > oo) and V(r) is any potential that 
can can be decomposed in a Fourier expansion. The method applies in general if V(r) includes various possible effects 
such as: lattice periodic potential, finite wave-length corrections, impurity scattering, etc; however as it will be lately 
argued, the impurity scattering is the most likely explanation for the recent experimental results. One important 
remark with relation to the impurity potential in Eq. Q is that to start with, it should only include the polarization 
effects produced by the combined effects of the Landau-Floquet states and the impurity potentials, the broadening 
effects produced by this potential are, as usually, included through the Kubo formula. Then we write V(r) as 

V(r) = V(r)-AV(r), AV(r) — V(r)W , (2) 

were W is the transformation that takes exactly into account the microwave-Landau dynamics, it is explicitly given 
in Eq. (|10|l . In the absence of the microwave radiation, W = 1 and V(r) vanishes. The impurity scattering potential 
V(r) is decomposed as 



(r-r,)}, (3) 



here r» is the position of the ith impurity and the explicit form of V(q) depends of the mechanism that applies under 
particular physical conditions, some examples will be consider in section (|IV|) . H{g w j is then written as 

H {B ^ = ^-U 2 , n = p + eA, (4) 

m* is the effective electron mass and the vector potential A includes the external magnetic and radiation fields (in 
the A — ► oo limit) contributions 



-r x B + Re 



— - exp{— iut} 

UJ 



(5) 



We first consider the exact solution of the microwave driven Landau problem, the impurity scattering effects are 
lately added perturbatively. This approximation is justified on the following conditions: (i): \V\/huj c << 1 and (ii): 
toTtr ~ uj c T tr » 1; r tr is the transport relaxation time that is estimated using its relation to the electron mobility 
/i = ertr/m*. As discussed in section (|IV|) both conditions are fully complied. 

The system posed by H{ B w j can be recast as a forced harmonic oscillator, a problem that was solved long time 
ago by Husimi 21 . Following the formalism developed in references22i2£, we introduce a canonical transformation to 
new variables P M ; /j = 0, 1, 2, according to 

Qo = t , Po = id t + ecf> + er ■ E, 

V^BQi = n y , VeBPi = n„ 

^f&BQi = p x + eA x + eBy , Ve~BP 2 = p y + eA y - eBx. (6) 

It is easily verified that the transformation is indeed canonical, the new variables obey the commutation rules: 
— [Qo, Po] — [Qij Pi] — [Q2, P2] — i] ah other commutators being zero. The inverse transformation gives 

x = l B (Qi -Pi), y = l B (Q2 - Pi) , (7) 
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where l B — y -^g is the magnetic length. The operators (Q2, P2) can be identified with the generators of the electric- 
magnetic translation symmetries'^. Final results are independent of the selected gauge. From the operators in Eq. 
we construct two pairs of harmonic oscillator-like ladder operators: (a\,a\), and (0,2, a\) with: 

ax = (Pi - iQi) , a 2 = -±= (P 2 - 1Q2) , (8) 

obeying: [ai,a{] = [02, a\] = 1, and [01,122] = [01,03] = 0- 

It is now possible to find a unitary transformation that exactly diagonalizates H{ Bu ,y, it yields 

W^H {BtU} W = u;J^ + a\ai\ =H , (9) 

with the cyclotron frequency uj c = eB/m* and the W(t) operator given by 

W(t) = exp{i?7 1 (3 1 }exp{^ 1 P 1 }exp{i772Q2}exp{^2P2}exp{-i J Cdt'} , (10) 

where the functions r]i(t) and £i(i) represent the solutions to the classical equations of motion that follow from the 
variation of the Lagrangian 

£ = y (vi + C1 2 ) + Cim + + d B [E x (Ci + m) + e v ( m + C 2 )] ■ (n) 

The explicit form of the solutions for r)i (t) and ^ (t) are given in the appendix IjVI A(l . 

Let us now now consider the complete Hamiltonian including the contribution from the V(r) potential. When the 
transformation induced by W(t) is applied the Schrodinger equation in becomes 

P *W = H ^ w > + V w (t)¥ W) , (12) 

where ¥ w > = W(t)V and 

V w (t) = W^Vi^W-^t) = Wi^V^W- 1 ^) - V{r) . (13) 

Notice that the impurity potential acquires a time dependence brought by the W{t) transformation. The problem is 
now solved in the interaction representation using first order time dependent perturbation theory. In the interaction 
representation = expliHot}^^^, and the Schrodinger equation becomes 

id t 9W={V w (t)} I ¥™ ) . (14) 

The equation is solved in terms of the evolution operator U(t), in such a way that ^^(t) = U(t — io)^j (*o)- The 
solution of the evolution operator in first order perturbation theory is given by the expression 

U{t) = l-if dt'\w^(t')V(r)W(t')\ , (15) 
J— 00 1 



that is explicitly evaluated in the appendix (|VI B|l . The interaction is adiabatically turned off as to — > —00, in which 

case the asymptotic state is selected as one of the Landau-Floquet eigenvalues of Hq, i.e. |\j>^ W ' ) (^o)) — * \l JL -k). The 
solution to the original Schrodinger equation in Eq. Q has been achieved by means os three successive transformations, 
which expressions have been explicitly obtained: 

|* M , fe (t))=W t exp{-iH t}U(t-t Q )\^k). (16) 

As discussed on the appendix (|VI A|) the index k labels the degeneracy of the Landau-Floquet states. Selecting the 
P-representation the dependence of the wave function on k becomes very simple, see Eq. I|47fl . for simplicity in what 
follows the index k will not be shown. The expression of the Kubo formula that will be derived in section (|IIIfl requires 
the knowledge of the matrix elements of the momentum operator II: 

(iJVl II; = (^\u\t-t ) [WniW^]jU(t-to)\v). (17) 
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Let us first consider the term inside the square brackets, using the explicit form of the operators in Eqs. ©. and 
(fTUfl it yields 



/eB(Qi -&), 



i = x, 
i = y. 



(18) 



If we now utilize the result for the evolution operator U given in the appendix (|VIBf) . we can explicitly work out the 
matrix elements of the momentum operator 



eB 



-iuicts; \ 



+ Ve~B~J2 e^+^-W A»(j). 



(19) 



Here the following definitions were introduced: = — £ v , a,j = bj = 1 if j = x and a,j = —bj 
A^(j) is given by 



A$(j)=*m* [pjSi,i+PjSi,-i] 



1 

71 



£„„ — uj c + uil — ir\ £ M „ + oj c + ujI — in 



-i if j = y, and 



(20) 



where q = ilsiqx — iqy)/V2j and the expressions for the functions pi, and „ are worked out as 



Pi 



el B E (~iuje x + u> c e y ) 
lo 2 - lu 2 + itoT rad 



92 



el B E (u c e x + iojty) 

LU 2 -U) 2 + iojT rad 



and 



= E / ^V( q )e-^D^ q) (A)' j t (| A | 



(21) 



(22) 



where Ji — Ji — <5/.o, J; being the Legendre polynomials, and D^ v (q) is given in terms of the generalized Laguerre 
polynomials in l|53|) . and 



A 



u) c l B eEo 



UJ (UJ 2 - UJ 2 + iujT ra d) 



[to (q x e x + q y 6y) + iuj c (q x e y - q y e x )j 



(23) 



It is important to notice that the subtracted term J; = J; — 5i t o, has its origin in the fact that the impurity potential 
Eq. (J2J includes only the dynamical effects, with the corresponding zero field-term conveniently subtracted. This 
procedure is justified because the broadening effects produced by V(r) are separately included via the Kubo formula, 
see appendix (I VI E|) . The subtraction J; = Ji — Sip becomes essential, otherwise the longitudinal resistance would be 
dominated by the / = term, producing incorrect results. 



III. KUBO FORMULA FOR FLOQUET STATES. 



In this section we shall develop the Kubo formula that applies when the dynamics includes Landau-Floquet states 
as those in Eq. i|ltj|) . We take the perturbing electric field to have the form E ext = Eq cos (fit) exp (— r]\t\). The 
static limit is obtained with £1 — > and n represents the rate at which the perturbation is turned on and off. The 
perturbing electric field is included in the vector potential, as we are interested in the linear response the perturbing 
potential has the form 



1 



n a p 



Eq 



sin (fit) exp (— rj\t\) 



(24) 



Besides the original Hamiltonian in JTJ , the complete Hamiltonian should include V ext and the part of the disorder 
potential (AV(r) — V(r)W) that was previously subtracted (see Eq. @). Hence the total Hamiltonian Ht is 
written as 



Ht = H + V e: 



AV(r) 



(25) 
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The disorder potential AV(r) will induce broadening effects, and it will be lately included. Then, the time-evolution 
for the density matrix p(t) obeys the von Neumann equation 

ih^- = [H + V exU p\. (26) 
ot 

Within the linear regime p is split in the sum p = po + Ap. The zero order term po must satisfy 

ih^ = [H,p ], (27) 

the conditions required to solve this equation will be established below. The first order deviation Ap then obeys 

dAp 



Of 



= [H,Ap} + [V ext ,p } . 



(28) 



We shall now apply to this equation the three transformation that were utilized in the previous section in order to 
solve the Schrodinger equation, hence in agreement with Eq. i|ltjfl . Ap is defined as 



Ap(t) = U\{t - t ) exp{iH t}W(t) Ap(t) W\t) exp{-iH t}Ui(t - <o) ■ 
In terms of the transformed density matrix Ap(t), Eq. I|28|) becomes 



ih 



dAp 
dt 



Vext i Po 



(29) 



(30) 



where V ex t and po are the external potential and quasi-equilibrium density matrix transformed in the same manner 
as Ap is transformed in Eq. I|29|l . The transformed quasi-equilibrium density matrix is assumed to have the form 
po = f{£n){lAi where /(£ M ) is the usual Fermi function and £ M the Landau-Floquet levels (see appendix I VI A() . 

It is straightforward to verify that this selection guarantee that the quasi-equilibrium condition in l|27|l is verified. The 
justification for selecting a Fermi-Dirac distribution in the quasi-energy states is presented in the appendix (|VI C|l . 
It is shown, that under experimental conditions (r w <C Tt r <C Ti n ), the elastic and inelastic relaxation processes can 
be neglected as compared to the external field effects. The solution of the Boltzmann equation yields, for a weak 
microwave intensity, a Fermi-Dirac distribution in the quasi-energy states. The expectation value of Eq. (|30|l in the 
\p) base can now be easily calculated using Eqs. ifTBI) . l(2*^|) and (|2*§1) . solving the resulting equation with the initial 
condition Ap{t) — > as t — > — oo, yields for t < 



( M | Ap(t) \u) = Ap(t) 



eE Q 



i(Q—ir])t' 



n(t') \v v ) + (o -f -n) 



(31) 



where the definition f^ u = f{£^) — f{£ v ) was used. Substituting the expectation value for the momentum operator 
given in (|f 9JI . the integral in the the previous equation is easily performed. The current density to first order in 



Tr 



A P (t)J{t) 



, the resulting expression represents 



the external electric field can now be calculated from (J(t,r)} 
the local density current. Here we are concerned with the macroscopic conductivity tensor that relates the spatially 
and time averaged current density j = (r w F)~ J Q " dt J d 2 x(J(t,r)} to the averaged electric field; here t w = 2tt/lu. 
The macroscopic conductivity can now be worked out, results for the dark and microwave-induced conductivities are 
quoted: 



n 



fl I' 



iAh 



aip5^ 



M.^+i 

£ tiv + ft — ir] £ u 



biv8 



IT) 



(32) 




t-f £ M „ + Lul + Q - ir\ J 

In these expressions the external electric field points along the x-axis. Hence, setting i = x or i = y the longitudinal 
and Hall conductivities can be selected. The denominators on the R.H.S. of the previous equations can be related to the 
advanced and retarded Green's functions G^(£) — 1/ (£ — £ lL ± irf). To make further progress the real and absorptive 
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parts of the Green's functions are separated taking the limit r/ — > and using lim v ^o l/(£ — irf) — Pl/£ + iir5(£), 
where P indicates the principal-value integral. As usual the real and imaginary parts contribute to the Hall and 
longitudinal conductivities respectively. In what follows details of the calculations are presented for the longitudinal 
microwave-induced conductivity, the corresponding dark conductivity expressions as well as the Hall microwave- 
induced conductance are quoted in the appendix (|VI D|) . Implementing the previous considerations and inserting a 5 
function, the longitudinal microwave-induced conductivity takes the form 

+ (n->-n)l 1 (34) 



Ah 



where ImG v {£) = ^ [G^(£) — G~(£)]. The static limit with respect to the external field is obtained taking f2 — > 0. 
In the case of the impurity assisted contribution an additional average over the impurity distribution has to be carried 
out, it is assumed that the impurities are no correlated, utilizing the explicit expressions for the velocity matrix 
elements in 1)21) |). the final result for the averaged microwave induced longitudinal conductance is worked out as 



^ J II.U 1 L 



]2 I d2 1 2 

r n mp i B i ^ 2 q y 



' J l (\A\)V(q)D l±v (q) 

(35) 

where ni mp is the two dimensional impurity density and the following function has been defined 

B®(s,e v ) = - 



-^{[f(£ + Ioj + £o) - f(£)} Im G v {£ + Ilo + £ Q )} 
d£ 



(36) 
£o=0 



The photoconductivity in (|35|l has a first contribution that depends on the p\ factor (independent of the impurity 
concentration), it represents the direct cyclotron resonance heating, arising when the W(t) transformation is applied 
to the momentum operator (see Eq. UlSj l). The impurity induced contribution (second row in l|35|ll takes into account 
the dynamics produced by the magnetic and microwave fields, combined with the resonant effect of the impurities; 
the information is contained in the complete wave function in Eq. 

The previous expression would present a singular behavior that is an artifact of the 77 — >■ limit. This problem 
is solved by including the disorder broadening effects. A simple phenomenological prescription is dictated by simply 
retaining a finite value of 77 that is related to the quasiparticles lifetime (77 = 27r/r s )2LS£. According to this prescription 
the density of states (DOS) of the /1 level would have a Lorentzian form ImG^E) = j^s^TrpjA ' ^ more formal 
procedure requires to calculate the broadening produced by, the so far neglected, part of the disorder potential AV(r) 
(see Eq. I|25|)'l: fortunately as explained in the appendix (|VI E() the calculation becomes equivalent to that carried out 
by Ando^ii and Gerhardts^i, so the density of states for the /z-Landau level can be represented by a Gaussian-type 

f orm 32 

ImG,{£) = y^exp [-(£ - £,f/{2Tl)] , r» = B^2. (37) 

The parameter /3 M in the level width takes into account the difference of the transport scattering time determining the 
mobility /x, from the single-particle lifetime. In the case of short-range scatterers r tr = t s and j3n = 1. An expression 
for fin , suitable for numerical evaluation, that applies for the long-range screened potential in 1)40(1 is given in appendix 
(|VI Efl . /3 M decreases for higher Landau levels. This property becomes essential to generate NRS, because they only 
appear for a narrow T^, a condition that is satisfied around the Fermi level in the case of large filling factors. 

Eqs. H35I36[) contain the main ingredients that explain the huge increase observed in the longitudinal conductance 
(and resistance), when the material is irradiated by microwaves. In the standard expression for the Kubo formula there 
are no Floquet replica contribution, hence ui can be set to zero in (I36II . if that is the case becomes proportional 
to the energy derivative of the Fermi distribution, that in the T — * limit becomes of the form S(£ — £f), and the 
conductivity is positive definite depending only on those states lying at the Fermi level. On the other hand, as a 
result of the periodic structure induced by the microwave radiation, B"' contains a second contribution proportional 
to the derivative of the density of states: -^ImG v {£ + lui). Due to the oscillatory structure of the density of states, 
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this extra contribution takes both positive and negative values. According to Eq. (|37|l this second term (as compared 
to the first one) is proportional to the electron mobility, hence for sufficiently high mobility the new contribution 
dominates leading to negative resistance states. The former observation becomes fundamental, because in agreement 
with experiment, our calculations show that NRS can only be induced in very high mobility samples (see Fig. yj. 

As it was mentioned in section {nj the present method applies in general if V(r) can be decomposed in its Fourier 
expansion e.g.: finite wave-length corrections, lattice periodic potential, impurity scattering, etc. The microwave 
radiation by itself only produces transitions between adjacent Landau levels (first term on the R.H.S of Eq. (|35|) ) 
leading to the cyclotron peak. In the case of a periodic potential the resulting spectrum will be dominated by the 
region q w 2ir/a, where a is the lattice parameter; for the experimental conditions Xb ^> a and the contribution is 
negligible. So we are led to analyze the impurity assisted mechanism as a plausible scenario to explain the strong 
oscillatory structure of the magnetoresistance. 



IV. RESULTS. 



The single-particle and transport relaxation rates induced by disorder are given b; 

m* f 2kF \V(q)\ 2 f 1 



,29 



1 f — '"imp 



dq 1 "-V xo X i q * , (38) 



r > Kh 3 k F Jo i/l - (q/2k F ) 2 ( 2k 



where V(q) is the Fourier transform of the impurity potential Remarkably, we have a consistent formalism in 
which: (i) the photoconductivity l|35|) . (ii) the relaxation rates l|38[) . and (iii) the level broadening Eqs. (|37|l and 1G0I) 
can all be consistently calculated once V(q) has been specified. 

For neutral impurities the potential can be represented by a short range delta interaction, the coefficient in Eq. @ 
corresponds to a constant that can be selected as V(q) = 2irh 2 a/m* . The expression in Eq. (|38|l is readily calculated 
to yield the same value for the single-particle and transport relaxation rates 

a n imp ' (39) 



(TV) (JV) ' * " "imp 

T tr ~s 

the upper index N labels the neutral impurity case. The evaluation of the photoconductivity l|35|l requires in general 
a time consuming numerical integration. However, for moderate values of the microwave radiation the transitions are 
dominated by single photon exchange, in the neutral impurity case a very precise analytical approximation can be 
explicitly worked out; see appendix (|VI F() . 

For charged impurities the Coulomb potential is long-range modified by the screening effects. Although electron 
motion is restricted to 2-dimensions, the electric field is three dimensional and there are contributions from the impu- 
rities localized within the doped layer of thickness d. The screened potential can then be represented in momentum 
coordinates by the expression^*! 

rr^ vtf e~ 9d e 2 m* 

V (fi = — i , g ■ 1TF = 7. Z2, 40 

m 1 + 777 2Tre e b h 2 

where the Thomas-Fermi (TF) approximation is implemented in order to calculate the dielectric function. Here 
ei, represents the relative permittivity of the surrounding media. The expression in (|4(J|) corresponds to a screened 
potential, that in real space has a r~ 3 decay for large r. The rates in Eq. (|38|l can be evaluated numerically, however 
an accurate analytical results is obtained observing that the decaying exponential in (|40|) causes the integral to die 
off for q » l/\d\, and the upper limit in the integral can therefore be set to infinity. Additionally for the relevant 
parameters (see below) the following conditions are observed k F » f/M and qrF >> 1/| rf| ; consequently, it is 
reasonable to drop the factor q/qTF in the denominator of l|4U|l and replace the square root in the denominator of Eq. 
(|38|l by unity. These simplifications yields for the transport relaxation rate 



1 irh 



_ (C) 

rf> 8m*(k F df nmp ' 



(41) 



the upper index C labels the charged impurity case. As expected, for charged remote impurities the single-particle 
lifetime differs from the transport lifetime, the approximated relation reads r t ^ ~ (2kpd) 2 Ts°^ . 

The parameter values have been selected corresponding to reported experiments^ in ultraclean GaAs / Al x Gai- x As 
samples: effective electron mass m* — 0.067 m e , relative permittivity ef, w 13.18, fermi energy £p = lOmeV, electron 
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FIG. 1: Longitudinal resistance, both total (continuos line) and dark (dotted line) as a function of the magnetic field. The 
figure also includes the Hall resistance (dashed line, p xy is rescaled by a factor 1/10). Results corresponds to neutral impurity 
scattering obtained with the approximated solution f appendix IVI Ffl and the selected parameters are: p « 0.25 x 10 7 cm, 2 /Vs, 
T ~ 1 K, f — 100 Ghz, \E\ w 2.5 V/cm, a 2 n\^ p = 5 x 10 6 cmT 2 . The values of the other parameters used in the calculations 
are discussed in the text. 



mobility /i Ri 0.1 — 2.5 x 10 7 cm 2 /Vs, electron density n = 3x 10 11 cm -2 , microwave frequencies / = 50 — 100 Ghz, 
magnetic fields in the range 0.05 — 0.4 Tesla and temperatures T ~ 0.5 — 2.5 K. The reported specimen is an 
5 mm x 5 mm square. Typical microwave power is 10 — 40 mW , however it is estimated^ that the microwave power that 
impinges on the sample surface is of the order 100 — 200 /jW, hence the microwave electric field intensity is estimated 
as \E\ w 1 — 3 V/cm. Using these values, it is verified that the weak-overlapping condition holds: uj c r tr ~ 100 — 1000. 

Recalling that /i = er tr /m*, one can use Eqs. 1|39[) an d l|41l) to determine the values of n imp corresponding to 
neutral or charged scatterers respectively. For example assuming fi w 2.5 x 10 6 cm 2 /Vs, one estimates for neutral 
scattering a 2 « 1 x 10 7 cm -2 . Although a and n\^ p are n °t separately fixed, one notices that the condition for 
weak disorder potential as compared to the Landau energy can be be expressed as V{q)/{l 2 ^hw c ) = 2ira <C 1; e.g if 
a ~ 0.01, then n^^ p ~ 10 11 cm -2 . For charged impurities and taking a value for the separation d between the impurity 
and the 2DES as d « 20 nm, yields nf/^ p ~ 1.5 x 10 11 cmT 2 . In this case the weak disorder condition takes the form 
V{q)/(l 2 B fvjj c ) ~ 7rexp(— 2-Kd/ls) <C 1, that is satisfied. A final remark is related to the radiative electron decay T ra d 
that determines the direct electron response to the microwave excitation (see Eqs. (|45|) and l|23() ). following reference 34 
Trod is related to the radiative decay width that is interpreted as coherent dipole re-radiation of electromagnetic waves 
by the oscillating 2D electrons excited by microwaves. Hence, it is given by V r ad = ne 2 / (6eom*c), using the values of 
n and m* given above it yields T ra d ~ 0.38 meV. 

Adding the dark and microwave induced conductivities, the total longitudinal a xx = + {(J xx ), and Hall a xy = 
o'xy + io'xy) conductivities are obtained. It should be pointed out that the interference between the dark and microwave 
contributions exactly cancels. The corresponding resistivities are obtained from the expression p xx = a xx / (a xx + tj 2 y ) 
and p xy = a xy / (<r xx + cr 2 y ). The relation a xy 3> o xx holds in general, hence it follows that p xx oc <r xx , and the 
longitudinal resistivity follows the same oscillation pattern as that of a xx . The plots of the total longitudinal and Hall 
resistivities as a function of the magnetic field intensity are displayed in Fig. ^ Whereas the Hall resistance presents 
the expected monotonous behavior, the longitudinal resistance shows a strong oscillatory behavior with distinctive 
NRS. The behavior of the complete p xx is contrasted with the dark contributions that presents only the expected 
Shubnikov-de-Hass oscillations. 

Fig. |3shows a comparison of the longitudinal resistivity as function ofui/u> c obtained for the case of neutral impurity 
scattering using both the approximated expression in the appendix (|VI Ffl as well as the result of the numerical 
integration Eq. (|35|l : the electron mobility is selected as pi = 0.25 x 10 7 cm 2 / V s . The approximated analytical 
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FIG. 2: Longitudinal resistance as a function of e = lo/lu c for neutral and charged impurities. Results for neutral impurities 
are obtained from the numerical integration (dotted line) and also using the analytical approximation discussed in appendix 
IVIFfl (dashed line) with the parameters: p ~ 0.25 x 10 7 cm 2 /Vs, a 2 — 5 x 10 6 cm -2 . The continuos line corresponds to 

the charged impurity case with parameters p m 2.5 x 10 7 cm 2 /Vs, n^ p — 1.5 x 10 11 cm -2 . The other parameters are the same 
as in figure Fig. Q 



result shows a good agreement with the one obtained from the numerical integration. It should be remarked that 
the approximated expression includes only one-photon exchange processes, while the numerical result includes the 
possibility of multi-photon exchange, hence it is concluded that for the selected electric field intensity (\E\ ~ 2.5 V/cm), 
the one-photon processes dominate. Results are also presented for the case of charged impurity scattering but for 
/i = 2.5 x 10 7 cto 2 /Vs. In spite of the very different nature of the two physical processes and that the mobility is 
increased by an order of magnitude in the charged case, it is observed that the results for the neutral and charged 
cases are very similar. The similarity of both results is based on: (i) the increase in the mobility is compensated by 
the factor /3 M in (|37|l giving a similar broadening value, (ii) although for the neutral case V(q) is constant over all the 
g-range of integration, whereas for charged case V(g) varies according to the expression in (|4U|) . the integral in <|35l) 
is dominated by the region in which q sa 2tt/Ib- 

One of the puzzling properties of the observed huge magnetoresistance oscillations is related to the fact that they 
appear only in samples with an electron mobility exceeding a threshold value pth- The phenomenon is absent in 
samples in which p is slightly reduced. This behavior is well reproduced by the present formalism. Fig. [3] displays 
the p xx vs. e = u>/uj c plot for neutral impurity scattering and three selected values of p. For p « 0.5 x 10 6 cm 2 /Vs 
the previously known, almost linear behavior p xx oc B is clearly depicted. As the electron mobility increases to 
p w 1.5 x 10 6 cm 2 /Vs, the magnetoresistance oscillations are clearly observed; however, NRS only appear when the 
mobility is increased to p ~ 2.5 x 10 cm /Vs. It is observed that p xx vanishes at e = j for j integer. The period 
and phase of the oscillations follow a pattern very similar to the one observed in experiments 2 ^, with minima at 
e = j + 5, and maxima at e = j — S, adjusted with S w 1/5. It should be pointed out that this value of 5 depends on 
the correct representation of the density of states, using a Lorentzian form instead of the Gaussian in (|37|l would give 
5 ~ l/10Si. Similar behavior is observed for the charged impurity scattering case, but with the mobility threshold 
increased approximately by and order of magnitude p t h ~ 2.5 x 10 7 cm 2 /Vs. The precise determination of \x t h depends 
of course on the selected values of the other parameters, mainly on the frequency and microwave intensity. 

The first reported experimentsiS*^ were carried out for a microwave radiation with transverse polarization with 
respect to the longitudinal current flow direction. However, it is clear that the presence of the magnetic field induces 
a chirality in the system, the model can be used to test these effects. Fig. 0] shows the results for different E— field 
polarization's with respect to the current. In Fig. 0Ji it is observed that the amplitudes of the resistivity oscillation are 
slightly bigger for transverse polarization as compared to longitudinal polarization. This result is in agreement with the 
recent experiment 4 , in which it is reported that the selection of longitudinal or transverse polarization produces small 
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FIG. 3: Longitudinal resistivity as function of e = uj/ui c for neutral impurity scattering and three values of the electron 
mobility: p w 0.5 x 10 6 cm 2 /Vs (dotted line), p w 1.5 x 10 6 cm 2 /Vs (dash-dotted line), and « 2.5 x 10 6 cm 2 /Vs (continuos 
line). In the two former cases the oscillations follow a pattern with minima at e = j + S, and maxima at e = j — delta, adjusted 
with 5 ~ 1/5. NRS only appear when p > ixtn ~ 1-5 x 10 6 cm 2 /Vs. The values of the other parameters are the same as in 
figure Fig. 



differences. However, we propose that the more significant signatures will be only observable for circular polarization. 
Selecting negative circular polarization (see Fig. QJa), the oscillation amplitudes get the maximum possible value. 
Instead, for positive circular polarization an important reduction of the amplitude is observed leading to the total 
disappearance of the NRS. These results are understood recalling that for negative circular polarization and u ~ lo c 
the electric field rotates in phase with respect to the electron cyclotron rotation. Based on the present results, it will 
be highly recommended to carry out experiment for circular polarization configurations. 

The present formalism can also be used in order to explore the non-linear regime in which multi-photon exchange 
plays an essential role. As the microwave radiation intensity is increased, the analytical approximation breaks down 
and the numerical expression in (|35|) with higher multipole (I) terms needs to be evaluated. In the explored regime 
convergent results are obtained including terms up to the I — 3 multipole. Fig. [SJi displays p xx vs. e plots for 
electric field intensities \E\ — 2.5V/cm and \E\ — 5V/cm respectively, the increase on the field intensity produces a 
corresponding increase in the minima and maxima of p xx , but apart of this, the qualitative behavior in both cases 
is similar. A further increase of the electric field intensity to \E\ = lOV/cm and \E\ = 30 V/cm (Fig. \5}p), takes us 
to the nonlinear regime in which a qualitatively new behavior is observed. For e > 2 the same NRS are observed, 
however the width of these regions increase to include practically all the range from e = j to e = j + \. Notice 
that the negative resistance minima has not a monotonous dependence on \E\, in fact for the strongest field intensity 
the minima approaches zero. Remarkably, for \E\ = 30 Vy cm and e < 2, new negative resistance states associated 
with transitions by two microwave photons are observed near e = | and e = |. The minima of these states are 
centered at e m i n = 0.52 and e m j„ = 1.52 respectively. Evidence of ZRS associated with multi-photon processes have 
been already observed by Zudov etal. (— ); they reported structures with maxima near e = 1/2 and e = 3/2, and the 
corresponding minima centered around e m in = 0.67 and e m in = 1-68 respectively. Dorozhkin 5 and Willett etal. ( 6 ) 
have also reported p xx minimum associated with e = 1/2. Although the exact position of the minima and maxima 
of p xx observed in Fig. are not localized at the same position reported by Zudov etal. (— ), the general pattern is 
very similar, supporting the interpretation as multi-photon processes. Clearly, a more systematic analysis and further 
experimental studies are necessary. 

Comparison with some other theoretical work is obliged. Previous work in referencesiSiiiiiSiil analyzed the effects 
of the microwave radiation on the electron scattering by impurities in the presence of a magnetic field. Durst et alH 
consider an out of equilibrium calculation, instead here a quasi-adiabatic approximation is implemented, assuming 
that the system is thermalized in those states characterized by the Landau-Floquet spectrum. The similarity between 
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FIG. 4: Longitudinal resistance p xx for neutral impurity scattering and various microwave E— field polarization's with respect 
to the current. In figure (a) the continuos and dotted lines correspond to linear transverse and longitudinal polarization's 
respectively. Figure (6) shows results for circular polarization's: left-hand (continuos line) and right-hand (dashed line). The 
values of the parameters are the same as in figure Fig. Q 



some results in the present work and those of Durst et aU&, suggest that departure from equilibrium is not significant 
for the studied phenomenon. The present formalism extends and explores the impurity assisted photoconductivity 
mechanism in detail. In this model the same disorder potential determines the broadening of the Landau levels, as well 
as the wave function that is used to evaluate the velocity matrix element. These matrix element are incorporated into 
a Kubo-like formula that takes into account the Floquet structure of the system. As previously mentioned, there are 
alternative models in which the leading contribution arises from the modification of the electron distribution function 
induced by the microwave radiation. According to Dorozhkin^ the negative resistance phenomena has its origin in 
a local population inversion that produces a change of sign of the (df /d£) term that appears in the conductivity. 
Although possible, the inversion of population requires rather strong microwave powers, which were not achieved in 
the experiments^. Indeed, the inversion population is expected to be produced when the microwave energy exceeds 
the Fermi energy (eE) 2 / (m* lj 2 ) > Ep (see appendix IVI C|) . clearly the estimated value for the threshold electric 
field E t h ~ 1000 V/ cm highly exceeds the experimental microwave fields E ~ 1 — bV/cm. An interesting alternative 
explanation based on the modifications that the microwave radiation produces in the distribution function was recently 
presented by Dmitriev et ulM^ and Kennett et al.et alm&. In these publications it is assumed that the inelastic-scattering 
processes give the dominant contribution to the collision term of the kinetic equation. However, as explained in the 
appendix (j VI (J[) under experimental conditions r w -C r tr <C r, n , and certainly the inelastic processes can be safely 
ignored as compared to the elastic processes. In fact we have presented an argument for a first approximation in 
which the distribution function is determined only by the microwave effects. It may be interesting for a future work, 
to add to the present formalism the effects that elastic processes produces to the distribution function. In any case, 
we consider that the present results taken together with those of referencesi2*ii*i2ii£ consolidate the explanation of 
the photoconductivity oscillations and negative resistance states in terms of the microwave-disorder mechanism. 



V. CONCLUSIONS. 



We have considered a model to describe the photoconductivity of a 2DEG subjected to a magnetic field. We 
presented a thoroughly discussion of the method that allow us to take into account the Landau and microwave 
contributions in a non-perturbative exact way, while the impurity scattering effects are treated perturbatively. The 
method exploits the symmetries of the problem: the exact solution of the Landau-microwave dynamics © is obtained 
in terms of the electric-magnetic generators (@ as well the solutions to the classical equations of motion l|llfl • The 
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FIG. 5: Nonlinear effects in the longitudinal resistance p xx for charged impurity scattering. Figure (a) includes the p xx . v.s. 
e — lu/ujc plots for electric filed intensities of \E\ = 2.5 V/cm (continuos line) and \E\ = 5 V/cm (dotted line). Figure (b) displays 
results for \E\ — 10V/cm (dotted line) and \E\ = 30 V/cm (continuos line). The parameter values are: p, « 2.5 x 10 7 cm 2 /Vs, 
TttlK, f = 100Ghz, n£2 = 1.5 x 10 11 cm~ 2 . 



spectrum and Floquet modes are explicitly worked out. In our model, the Landau-Floquet states act coherently with 
respect to the oscillating field of the impurities, that in turn induces transitions between these levels. Based on this 
formalism, a Kubo-like formula is provided, it takes into account the oscillatory Floquet structure of the problem. 
It should be stressed that the disorder potential is conveniently split (see Eqs. J2J), i|25|) ) in such a way that it 
contributes both to the matrix elements of the velocity operator, as well as to the broadening of the Landau levels. 
Hence, we have a consistent formalism in which: (i) the photoconductivity (|35[) . (ii) the relaxation rates l|38|l . and 
(in) the level broadening Eqs. (j37(l and (|60[) can all be consistently calculated once the disorder potential has been 
specified. 

The expression for the longitudinal photoconductivity Eq. (1351) contains the main ingredients that explain the 
huge increase observed in the experiments. As explained in section l|IIIj) . the standard expression for the Kubo 
formula at low temperature is dominated by the states near to the Fermi level. On the other hand, as a result of the 
periodic structure induced by the microwave radiation the term i?W contains a second contribution proportional to 
the derivative of the density of states: jglm G U (E + Iuj). Due to the oscillatory structure of the density of states this 
extra contribution takes both positive and negative values. According to Eq. (|37|l this second term is proportional to 
the electron mobility, hence for sufficiently high mobility the new contribution dominates leading to negative resistance 
states. This allows us to explain one of the puzzling properties of the observed huge magnetoresistance oscillations, 
related to the fact that they appear only in samples with an electron mobility exceeding a threshold. This result is well 
reproduced by the present model, for the selected parameters, NRS emerge when the condition p, > 2.5 x 10 6 cm 2 /Vs 
(short-range disorder), and p, > 2.5 x 10 7 cm 2 /Vs (long-range disorder) are satisfied. The oscillations follow a pattern 
with minima at e = j + 5, and maxima at e = j — S, adjusted with S sa 1/5. These results are in reasonable good 
agreement with the observation of Mani etalm^, they reported a similar pattern with 5 ~ 1/4. 

An interesting prediction of the present model is related to polarization effects that could be possible observed 
in future experiments. While the results for the cases of linear transverse or longitudinal polarization's show small 
differences, the selection of circular polarized radiation leads to significant signatures. The maximum possible value 
for the oscillation amplitudes of p xx appears for negative circular polarization. Instead, positive circular polarization 
yields an important reduction on the oscillation amplitudes and the total disappearance of the NRS. This result can 
be understood, if one recalls that for negative circular polarization and u> ~ uj c the electric field rotates in phase with 
respect to the electron cyclotron rotation. The present results call for the importance to carry out experiment with 
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circular polarization configurations. 

An analysis was presented in order to explore the non-linear regime in which multi-photon exchange plays an 
essential role. The results suggest the existence of new NRS (which are expected to develop into ZRS) near e = 1/2 
and e = 3/2, these states correspond to two-photon exchange processes and are in reasonable agreement with the 
reported experimental results. 

Some final remarks are related to the limitations and possible extensions of the present work. In a first approximation 
we have not included the contribution of the elastic processes to the kinetic equation that determines the electron 
distribution. However, it will be interesting to extend the present calculations to include, not only the dynamical 
effects produced by the impurity on the electron wave function, but also the modifications that they produce in the 
distribution function. 



VI. APPENDICES. 



A. Microwave driven Landau problem. 



Eq. © defines a canonical transformation from the variables {t, x, y',Po,Px,Py} to {Qo, Qi, Qi\ Po, Pi) P2}, in terms 
of the new variables the Schrodinger equation Q (without impurity potential) takes the form 



Pn* = 



huj, 



Qi 



el B E x (Qi - P 2 ) + el B E y (Q 2 - P x ) 



(42) 



The action of the transformation W defined in Eq. (|10fl over the (Q^, P M ) variables can be easily calculated as: 



WQ W^ 
WP W* 
WQiW^ 
WQ 2 W^ 



Qo, 

Q1 + C1. 
Q2 + C2, 



C1P1 + V2Q2 + C2P2 - Cuyi - C2V2 + C, 

WP X W^ = Pi — 77 l5 

wp 2 w^ = p 2 - % . 



(43) 



It can be verified that when the W transformation is applied to Eq. I)42|l . the second and third terms in the right 
hand side exactly cancel with all the terms that appear in the expression for WP W^ (except P ) if £ is identified 
with the Lagrangian in Eq. (|1 1|) an the functions r\i and Q satisfy the differential equations that follows from the 
variation of C: 



i}2 



elsEy, 



Ci 

6 



u c rji = -ekEy, 
■■ -elsEx . 



(44) 



But, these are exactly the classical equations of motion that follow when the variational principle is applied to the 
Lagrangian in Eq. I|ll|) . Hence the W operator transforms the Schrodinger equation (|42fl to the Landau eigenvalue 
problem with the Hamiltonian given in Eq. @. 

For the electric field consider the expression in Eq. (JSJ. It is then straightforward to obtain the solutions to Eqs. 
B44J1 . adding a damping term that takes into account the radiative decay of the quasiparticle, they read 



Vi 

Ci 



cIbEq Re 
cIbEq Re 



-iue x + uj c e y - t 



rad 



uj c e x + ibJty jut 



- iujT. 



rad 



V2 
C2 



cIbEq Re 
—elsEo Re 



e v e" 



(45) 



According to the Floquet theorem the wave function can be written as \I/(i) = exp (— iS^t) <j>n(t), where 4>^(t) is 
periodic in time, i.e. 4>^(t + tJ) = 4>^{t). From Eq. (|1U|> it is noticed that the transformed wave function ^> w = 

contains the phase factor exp (if Cdt' J . It then follows that the quasienergies and the Floquet modes can be deduced 



if we add and subtract to this exponential a term of the form £ JJ* Cdt' . Hence, the quasienergies can be readily read 
off 



f(0) 



P(0) 



huj r 



£rad 



e 2 El [1 + 2m c Re(e* x e y )/u 
2m* Uu-ujcf + r^ 



(46) 
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here Sjf are the usual Landau energies, and the induced Floquet energy shift is given by the microwave energy £ ra d- 
The corresponding time-periodic Floquet modes in the (Pi, P 2 ) representation are given by 

* M ,*(P) = exp{-ism(2cut) F(lu)}^(P 1 )5{P 2 - k) , (47) 

here ^ (Pi) is the harmonic oscillator function in the Pi representation 

1 



(Pi) = (Pi | M ) = 



3-^/ 2 PT m (Pi) 



and (Pi) is the Hermite polynomial and the function F(u>) is given as 



eEnl 



2 r 



^ - ^ + 2u^ x - 2u^ + (2^ - u,W c + .j) 



(48) 



(49) 



The wave function 147|) depends on the Landau (/1), and center guide (k) indexes; however the spectrum (|4(j|| is 
degenerate wit respect to k. It is important to notice that F(u>) appear in the wave function phase, that depends only 
on time, hence its contribution to the expectation value of the momentum operator cancels exactly. Thus, contrary 
to what it is claimed in reference^, the effect of the Floquet dynamics (without including an extra effect such as 
impurity scattering) can not account for the explanation of the ZRS observed in recent experiments. 

B. Impurity induced transitions. 

In this appendix we consider the first order solution of the evolution operator U (t) given by 



U(t) 



1 



dt' 



W\t')V(r)W(t') 



(50) 



The effect of the transformation induced by the W operator over the impurity potential can be easily evaluated 
considering the effect over the Fourier decomposition of V(r) given in Eq. Recalling that the x and y coordinates 
are written in terms of the new variables (Qi, Pi, Q 2 , P2) by means of Eq. (JJJ, and utilizing the transformation 
properties of the (Qi,Pi) operators in (|43|l . it is readily obtained 



W\t) exp{iq ■ r}W(t) = exp{il B (q x P 2 - q y Q 2 )} x 

exp{-i/ B (q x Qi - q y P\)} x exp{iZ B [q x (Ci + m) + Qy (C2 + m)]}- 



(51) 



Using Eqs. (|45H the third exponential in the previous equation can be recast in a compact form as 
exp{— iRe (A exp(iu)t))}, with A given in Eq. (|23p; this expression can be expanded aa^ 



exp{— iRe (A exp(iujt))} = 



ilAI 



Jim), 



(52) 



with Ji the Legendre polynomials. For the the second exponential notice that once that Q\ and Pi are replaced 
by the raising and lowering operators given in Eq. ijHJl, one is lead to evaluate the matrix elements of the operator 

D (q) = exp (qa\ — q*ciij that generates coherent Landau states. A calculation yields 



{q) = (v\D{q)\n) = 



(53) 



where are the generalized Laguerre polynomial. With all these provisions the matrix element of the solution of 
the evolution operator in (|50|) can be worked out as 



(n\U{t)\v)=8 llv -Y J 



— c® 

C 1 7 UK' 

tfj,v ~ IT] p 



(54) 



the explicit expression for C^l was given in 11' 21) 
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Microwave-driven distribution function. 



Within the time relaxation approximation the Boltzmann equation can be written as 



df 
dt 



dp 



(eE + ev x B) 



f-fr f-fr 



where fp is the Fermi-Dirac distribution and we distinguish between the elastic rate t 4 ~ and inelastic or energy 
relaxation rate . Under experimental conditions: t w <C r tr <C Ti n , and certainly the inelastic processes can be 
safely ignored. Furthermore, due to the ac-electric field Q, the L.H. S. of the previous equation is estimated to be 
of order f /t u \ hence, in a first approximation the elastic scattering contribution can also be neglected. The resulting 

Vlasov equation has the exact solution f(p, t) = fp (p — m*v(t)), where the velocity v(t) = £ij solves exactly 

the same classical equations of motion as given in (|44f) , and the initial condition is selected as / — * fp as the external 
electric field is switched-off. In particular it is verified that m*\v(t)\ 2 /2 = £ ra d coincides with the Floquet energy shift 
produced by the microwave radiation l|4tj[l . The steady-state distribution, evaluated at the Landau energy £ = £^ 
obtained by averaging fp (p — m*v(t)) over the oscillatory period 



(0) 



(If) = - 



(I 



fp ( 4°) + £ rad + 2 cos uj c t J 4 0) £ r 



For the experimental conditions it is verified that £ ra <i <C , thus expanding to first order one finds (fp) rj 

fp + Sradj = /f(^)- Hence, it is verified that a rapid relaxation of the Fermi distribution to the quasi-energy 

states is a reasonable assumption. The arguments presented in this appendix have been introduced by Mikhailo\«2* in 
order to explore the possibility that the microwave radiation leads to a population inversion; however, it is concluded 
that it would require a rather high microwave intensity £ ra d > £f- 



dt. 



D. Dark and Hall conductivities. 



In section (|III|) it was explained in detail the method to obtain the final expression for the microwave induced 
magnetoresistance Eq. (|35|) . Working along a similar procedure the expression for the remaining conductivities are 
worked from equations and First we quote the longitudinal dark conductance 



l=^I> / d£ImG^ (£) ^IrnG^ {£ + uj c ) , 



— 



(55) 



whereas the dark Hall conductance is given by 



M / d£ Im G p {£) [f {£» - u c ) - f (£)} P 



(£-£, 



(56) 



where P indicates the principal- value integral. The impurity assisted contributions require an additional average over 
the impurity distribution, it is assumed that the impurities are no correlated. The final result for the microwave 
assisted longitudinal conductivity was quoted in Eq. I|35|) . Following a similar procedure the microwave assisted Hal 
conductivity is calculated to give 



¥ / df E E /m ^ ( £ ) if ^ - f ft*)] (<w (pipi 



1 / nip 



d\ T(q) 



J l (\A\)V(q)D liu (q) 



were the function T(q) is defined as 

T(g)=c^| 



+ p\p25i-i) 



(E + wl-E,,) \{£ + - £„)- tu 2 c \ 2 



(57) 



(58) 
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E. Landau density of states. 



A detailed calculation of the density of states incorporating all the elements that contribute to the system under 
study is beyond the scope of the present paper. However, it can be argued that the expression given in Eq. (|37|) for 
the DOS is expected to be a reasonable selection under some consistent approximations. Le us consider the Green's 
function associated with the Hamiltonian H + AV, where H is given in and Ay = W^VW is the subtracted 
part of the disorder potential l|25[) . As explained in section (|lll|l . the Kubo formula (|32I33|) was deduced using the 
wave function obtained after the three transformations in Eq. (|16Jl are applied to the Landau states. Hence, AV 
is transformed according to A~V(t) = ll\(t - t Q ) e-np{iH t}W(t) AV W\t) evp{—iH Q t}Ui(t - t ) (see Eq. <P|) 1. 
Notice that: (i) the W transformation cancels exactly (ii) both AV and the first order correction to Ui(t — to) are 
proportional to V, hence considering linear terms on V we can set Ui(t — to) ~ lj finally when evaluated in the 
|/Lt) base and neglecting inter-Landau mixing, the contributions from exp{— iHot} cancels out. Hence AV rs V{ r )i 
and the problem under consideration reduces to evaluate the density of states produced by a magnetic field and a 
disorder potential V(r) of the form given in Eq. ©; but this is precisely the problem considered some ago time by 
by Ando^ . and Gerhardts^I; the density of states is well represented by the Gaussian expression in (|37[) . and the level 
broadening neglecting couplings between different Landau levels is taken from reference^: 

Tl=8*nlln mp J ^1 f^ V (r)V(r')[D^(\r-r'\(V2l B ))]\ (59) 

where is given in (|53|) . For the delta short range scatterers the previous expression is readely evaluated yielding 
the result in l|37|) with [3^ = 1. In the case of the charged impurity disorder, after the substitution of the Fourier 
decomposition © and using Eq. (|40|) it is verified that is again given by the expression in Eq. 1)3 7|) , but the factor 
is given by 



= 16tt (k F df I q ^ [D^ {q)Y dq . (60) 



exp 



2 



1 



V2q 
Ib Qtf 



Previous analysis of the Landau level broadening were carried out, for example^, for a Gaussian potential V(r) ~ 
e r l d , however as mentioned in section (|IV|) the actual situation corresponds to a screened potential that in real 
space has a r~ 3 decay for large r. As mentioned in section ill I II) the value of /3 M decreases for higher Landau levels. 
For example for the selected parameter, we have: fto = 108, /?3o = 14, (3^ = 11. 



F. Approximated one-photon exchange photoconductivity. 

The microwave induced longitudinal (|35|l conductivity requires the numerical evaluation of a time consuming integral 
given by 



Si = J d 2 qK(q) 



2 

J l (\A\)V(q)D^(q) . 



However, if we consider the regime of moderate microwave intensity and assuming neutral impurity scattering, a very 
useful analytical approximation can be worked out. For neutral impurity scatterers, the potential is assumed to be 
of the short range delta form, hence the Fourier coefficient in © is given by V(q) = 2irh 2 a/m* . The D^(q) term 
contains an exponential factor that represents a cut off for large q. Then according to Eq. (|23(l for moderate values 
of the microwave electric field the A term is small and the leading contributions arises from the I = ±1 factors that 
correspond to the single photon exchange contribution. Using the approximation J\{z) sa z/2 one is lead to evaluate 

S±1 = ^ f d 2 qK{q)\A\ 2 \D llv {q)\ 2 . (61) 

TO* ./ 



The angular integration is straightforward, while the integral over the q = \ q 2 + q 2 , leads, after a change of variable 



v 

„2 



£ = q to an integral of the form 

/>oo 

d £ e -i^-v+2 (L^ U ) 2 , (62) 
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that is explicitly evaluated with the help of the recurrence relation xV^ = (2n + k + 1)£„ — (n + k)L k l _ 1 — (n + l)L k l+1 
and the integral^ 



/>oo 

/ d^ k L k n L 
Jo 



k rk rk _ (Tl + k)\ 



n i 



Smn- (63) 



The final result reads 



where 



5 ±i = "f.T 2 '^- [^(1 + 2| e ,| 2 ) + c c 2 (l + 2| £:c | 2 ) - 8^ c Im (e* x e y )] , (64) 



V=6l/x+^ ) U+iT) +(M-^ + i (65) 



These expressions greatly simplifies the numerical calculations, and as it is discussed in section l|IV|) they provide a 
very accurate approximation to the exact result. 
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